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SUMMARY 


The induced drag polar is developed for wings capable of attaining 
extremely high loadings while possessing an elliptical distribution of 
circulation. This development is accomplished through a theoretical 
investigation of the vortex-wake deformation process and the deduction 
of the airfoil forces from the impulse and kinetic energy contents of 
the ultimate wake form. The investigation shows that the induced 
velocities of the wake limit the maximum lift coefficient to a value 
of 1.9^ times the wing aspect ratio, for aspect ratios equal to or 
less than 6.5, and that the section properties of the airfoil limit the 
lift coefficient to 12.6 for aspect ratios greater than 6.5- Relations 
are developed for the rate of deformation of the vortex wake. It is 
also shown that linear wing theory is applicable up to lift coeffi- 
cients equal to 1.1 times the aspect ratio. 

INTRODUCTION 


The development in recent years of advanced methods of boundary- 
layer control has made possible the attainment of very high lift coeffi- 
cients with two-dimensional airfoil sections. Use of boundary- layer 
control on three-dimensional wings has likewise resulted in large 
increases in the maximum sectional circulations by preventing flow 
separation at high angles of attack. For a wing of finite span, how- 
ever, this increased circulation is necessarily accompanied by an 
increasingly powerful trailing vortex system. The velocities induced 
in the vicinity of the wing by this system can become quite large, 
with a pronounced effect on the wing aerodynamic force, especially on 
the lift component. Successive increases in circulation produce 
diminishing increases in wing lift. 


1 The basic information presented herein was a part of a thesis 
entitled "The Limit of Circulation Lift on Airfoils of Finite Aspect 
Ratio" which was offered in partial fulfillment of the requirements 
for the degree of Master of Aeronautical Engineering, University of 
Virginia, Charlottesville, Virginia, May i960. However, some addi- 
tional material on the effects of wake deformation on longitudinal 
stability has been added. 


2 


From a design and performance standpoint it thus becomes of con- 
siderable importance to determine these wake effects in a quantitative 
manner. It is customary practice in conventional wing design to employ 
linear wing theory; this is justified since normal wing stall limits 
the maximum lift coefficient to relatively low values. For wings 
intended to operate at very high lift coefficients, however, the ques- 
tion arises as to the limit of the valid use of linear theory, since 
the linearizing assumptions are no longer justified for large deforma- 
tions and displacements of the wake. 

The present investigation has as its object, therefore, the quanti- 
tative estimation of the effects of the vortex wake deformation on the 
wing aerodynamic forces, with view to establishing the limiting lift 
coefficient for linear- theory application, and also the nature and magni- 
tude of the deviation from the linear prediction at the higher lift 
coefficients. This development is accomplished by considering in some 
detail the rolling-up process of the wake and by deducing the aero- 
dynamic forces from the impulse and kinetic energy of the final wake 
form. The results are presented in the form of the wing- induced drag 
polar, for all aspect ratios. As a corollary, for use in stability 
analyses, the rate of rollup of the vortex sheet is investigated and 
an expression is developed for the distance required for complete 
rollup, for a given lift coefficient and aspect ratio. 

The basic analysis is carried out for the specific case of a con- 
ventional wing which maintains an elliptical distribution of circula- 
tion at all lift coefficients. The various wing- propulsion systems - 
such as the jet- flap wing - are excluded. The elliptical distribution 
is chosen because it yields results which are generally applicable to 
the more common planforms where the larger part of the vorticity is 
shed near the wing tips. The methods of the analysis are, however, 
general and may be applied to other distributions, although the calcu- 
lations may be more involved. 

Throughout the analysis the air is treated as an inviscid, incom- 
pressible fluid. Inasmuch as the attainment of high lift coefficients 
by actual wings requires the removal or reenergization of the boundary- 
layer air, the production of a large viscous wake by the airfoil is 
precluded and viscous effects are negligible. In the sequel, therefore, 
the wake referred to is the vortex wake associated with the wing circu- 
lation gradient and not a viscous wake. 



SYMBOLS 


aspect ratio. 


b£ 

^w 


cross-sectional area of the ovoid wake 

circle radius (Joukowski transformation); ellipse major 
semi axis 


wing span; ellipse minor semiaxis 

D, 


induced drag coefficient, 


4 * 


lift coefficient, 



wing chord 
induced drag force 


vortex sheet rollup distance (along X-axis) 

vector element of surface area 

vector element of arc length 

total wake kinetic energy 

kinetic energy per unit length of subcore 

wing aerodynamic force vector 
— > 

component of F due to momentum changes 
resultant force on bound vortex system 


impulse per unit length of wake 
total wake impulse vector 



J 

L 

l 

M' 

n 

P 

— > 
q. 

R 2> r 


vortex moment of inertia about system center of gravity 
lift force 
generalized length 

hydrodynamic mass per unit length of wake 
number of turns per unit length of vortex helix 
denotes specific point in space 
local velocity vector of induced field 

velocity of center of gravity of generalized vortex system 
radial distances from vortex core centers 


r 


o 


r l^ r 2 


S 

s w 

t 

U 


u 


— > 


V,v 


Vr 


radius of fully developed vortex core 


R 1 

normalized radial distances from vortex centers, — and 

b 

respectively 


general surface area 


wing area 


time 

general stream velocity 
x-component of induced velocity 
free- stream velocity, (vector) 
resultant velocity at wing section 


v y-component of induced velocity 

Vrp circumferential velocity 

W complex potential, 0 + i\|/ 

w z- component of induced velocity 
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w 

— » 
X 


y 

z' 

z 

-» 

z 

a 


cm 


r 



2 

L 


0 

* 

-* 


velocity of center of gravity of vortex sheet 
unit vector along positive X-axis 

lateral coordinate of center of gravity of vortex system 
complex variable in Y-Z plane, y + iz 

vertical coordinate of center of gravity of vortex system 
unit vector along negative Z-axis 
geometrical angle of attack 
absolute angle of attack 


circulation, 



q.ds 


circulation around median section of wing 

inclination angle of final wake 
dummy variable of integration 
generalized wake inclination angle 

r D 

vortex strength, — 

2rt 

fluid mass density 


dynamic pressure 
the wake vortex sheet 

time required for complete rollup of vortex sheet 

velocity potential 

eccentric angle, cos -1 — 

b 

stream function 

surface vorticity of vortex sheet 
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Subscripts: 

c property of vortex core 

i index number 

S property associated with wing surface 

x,y, z components along X-, Y-, and Z-axis, respectively 

L property associated with vortex sheet 

Superscript: 

P reference point at which velocity is computed 

THE WAKE- DEFORMATION PROCESS 
Wake- Induction Effects 

Induced- velocity effects .- Consider a wing in steady flight 
operating with an elliptical distribution of circulation 


r(y) = rJii - 


( 1 ) 


Emanating from the trailing edge is a vortex sheet whose intensity 
immediately behind the trailing edge is 


dT 

d-y 


= -r. 


ir 


( 2 ) 


2yV 

b / 


This sheet may be considered a superposition of an infinite number of 

n P 

vortex filaments each of strength — dy as shown in figure 1. 

dy 

At any point P in the flow field the total velocity (relative to 

the fixed airfoil) is the sum of the free- stream velocity V, the 
-> P 

induced velocity q^ associated with the bound vorticity, and the 

P 

induced velocity q^ due to the wake vorticity. Vectorially, 
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L ' 
1 
1 

5 

6 


•-W 


-> P 
1 


— > 

= V p 


_1_ 

4rr 




( 3 ) 


This velocity field also applies to the vortex sheet and since the 
field is not, i n general, uniform, the vortex sheet undergoes deforma- 
tion as it passes downstream. Since the vorticity distribution of the 
wake determines the effective velocity field in the vicinity of the 
airfoil, by equation (3), this wake deformation is reflected as a change 
in the aerodynamic force vector acting on the wing. As is evident from 
equations (2) and ( 3 ), a large- span wing with small F Q (small angle 

of attack) will produce a relatively weak wake and the rate of deforma- 
tion will be small; thus, the wake will maintain its sheet form for a 
considerable distance downstream. 


High circulation effects .- If now the elliptical circulation loading 
form as given” by equation (l) is maintained while T Q is increased, the 

wake will increase in intensity according to the Helmholtz law of vortex 
continuity. The large induced velocities associated with the high inten- 
sity cause a rapid deformation of the vortex sheet which, in addition 
to being displaced downward, begins to roll up from the edges. The down- 
ward progression of the wake causes it to be inclined to the free- stream- 
velocity direction by an angle 0 = 9(x). The rollup process causes the 
vorticity of the initial sheet to become increasingly concentrated into 
two vortex cores, and the resulting alteration of the induced field 
decreases the downward progression of the wake so that 0, in general, 
decreases in the downstream direction (d0/dx negative). A generalized 
sketch of this variation of the wake angle is shown in figure 2. Here 
the angle 0 is the inclination of the forming vortex cores. Far 
downstream, where the wake has reached a stable state of deformation 
(has become fully rolled up) the inclination angle is denoted by 
5 = e(oo). 

The wake deformation has two primary effects on the flow at the 
wing. First, the result of the inclination of the wake is that the 
P 

velocity q induced at the airfoil is rotated forward as it increases 

in magnitude (fig. 3)- Second, the concentration of the vorticity into 
cores causes a variation of the magnitude of the downwash across the 
span. Since the wake intensity increases with r o so does the rate at 

which the sheet rolls up, and also, the angle of inclination of the wake 
increases. Thus, the downwash variation across the span becomes a func- 
tion of r 0 , for a wing of given span. For a wing with given r o , as 

the span is decreased the intensity of the wake increases and addition- 
ally the wake becomes more narrow, so that the deformation will proceed 
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more rapidly for a low- aspect-ratio wing than for one with high aspect 
ratio. Thus, the form of the downwash distribution is also a function 
of aspect ratio. 

The inclination of the wake introduces at the wing an induced- 

— > 

velocity component opposing the free- stream velocity V; thus, the 

— ) 

resultant effective velocity vector V R existing at a section of the 

wing decreases in magnitude and is rotated downward as the section cir- 
culation is increased. Both of these effects lessen the increment in 
lift which accompanies the circulation increase. *L 

In order to obtain an increase in circulation it is necessary to 
increase the absolute angle of attack of the section with respect 

— ► 

to Vr. However, as the geometrical angle of attack of a section is 
increased the velocity decreases in magnitude and rotates in a 

direction so as to reduce the absolute angle of attack; therefore, the 
variation of T with geometrical angle of attack becomes nonlinear in 
the high circulation range so that the slope dT/da decreases. Also, 

since the resultant velocity varies across the span, = Vp(y), due 

to the wake deformation, and the form of this velocity variation varies 
with r Q , the nature of the variation of sectional lift with angle of 

attack becomes quite complex. 

In the foregoing discussion the circulation distribution P(y) 
has been specified as being elliptical. In order to maintain such a 
distribution physically, it would be necessary to carefully tailor the 
wing geometry so as to accommodate the resultant- velocity variation 

V R (y). Since V R (y) also varies with T Q , it is clear that the wing 

geometry must be variable if the elliptical circulation distribution 
is to be maintained. Conversely, if the wing geometry is fixed, the 
form of r(y) must change with r Q . Of course, from this simple 
analysis no quantitative statements can be made about the nature or 
rapidity of the resultant- velocity variation. 

— ^ 

For the elliptical circulation loading, the value of Vp is uni- 
form across the span of an elliptical wing for small values of r Q . 

However, for large values of r o or small aspect ratios the rapid 
rollup of the vortex sheet will tend to change the Vp distribution. 

Hence even though an elliptical distribution of circulation may exist 

— > 

on a wing the lift distribution, which is dependent on V R , may there- 
fore be quite different in form. 


i H i— I Lf\VO 
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The results of this simple qualitative analysis indicate that the 
wake effects will set a limit on the maximum lift coefficient of a 
finite- span wing. In order to obtain a quantitative estimate of these 
effects it is necessary to consider the details of the wake deformation. 


Development of the Final Wake Form 

Consider now a straight wing of any aspect ratio operating with 
the elliptical circulation loading of equation (l). It is prescribed 
L that the loading shall remain elliptical at all values of lift coeffi- 

cient. This condition may be physically realized by suitably varying 
the sectional camber of the wing for each value of Cl or by properly 

twisting the wing. The intensity of the trailing vortex sheet at the 
trailing edge is given by equation ( 2 ). The flat vortex sheet is an 
unstable form of vorticity and, for reasons to be discussed later, pro- 
ceeds to roll in from the edges and to travel downward in the center so 
that relatively soon after its formation it appears somewhat as shown 
in figure 4. This rolling- up and displacement process continues until 
some distance downstream the original vortex sheet has separated in the 
middle and become coiled up into two cores of concentrated vorticity 
which are inclined by an angle 5 to the free- stream flow. These 
approximately circular cores extend on downstream to infinity where 
they join the ends of the starting vortex. 

If now the form and properties of the final wake can be quantita- 
tively determined, the aerodynamic force on the wing can be determined 
by suitable application of the general energy and momentum theorems. 

In order to establish the final wake form the system of reference axes 
shown in figure 5 is used, where the origin is located at the midspan 
of the trailing edge and the X-axis is parallel to the cores of the 
final wake, which is inclined to the free- stream velocity by the 
angle 5. The axes move with the airfoil so that the flow is steady. 
The positive directions of the induced- velocity components u, v, 
and w are as indicated. (The components of V are V cos 5 and 

V sin 5 in the x- and z-direction, respectively (fig. 6).) 

Under the action of the total velocity field the vortex filaments 
comprising the initial vortex sheet assume a helical form as they wind 
about the growing cores (fig. 4). Far downstream where the rolling up 
may be considered complete, this system of helical filaments approxi- 
mates what might be called a pair of vortex solenoids. These fluid 
solenoids possess the same induction properties as their electrical 
counterparts, with proper accounting, of course, for the fluidity of 
the cores. Far downstream the flow inside the cores has the components 

V cos 5 + u, v, and w - V sin 6, while outside the cores it has the 
velocity components V cos 6, v, and w - V sin 5 in the x-, y-, and 
z-direction, respectively. 
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Considering now a narrow strip of the flat vortex sheet taken near 
the trailing edge of the wing (as in fig. 4), if the rate of deformation 
of the sheet in the longitudinal direction is small, the wake may he 
established to a good approximation by considering the rolling up as a 
two-dimensional process and tracing the deformation of the two- 
dimensional strip as it moves downstream. At the higher deformation 
rates the assumption of a two-dimensional rollup appears less accurate 
since the shape of the wake on either side of the strip will vary some- 
what from that of the strip. However, the assumption is made that up 
to the maximum lift of the airfoil the rate of deformation is suffi- 
ciently low so that the two-dimensional deformation may be used to con- 
struct a first-order approximation of the wake form. A justification 
of this assumption lies in the fact that the primary concern is not 
with an accurate description of the entire sheet deformation but with 
only the relatively small region where the sheet is adding to the core. 
With this assumption the cross section of the wake at any distance 
downstream is that of the deformed two-dimensional strip. 

Properties of two-dimensional vortex systems .- The following two 
theorems which will be of subsequent use are derived in detail in 
reference 1: 

Theorem 1: If a system of two-dimensional vortices exists on one 

side of a flat bounding wall the center of gravity of the vortex system 
will translate parallel to the wall with the velocity 


<*0 


P iqi F R 

Pi 'Ppi 


( 4 ) 


where lb indicates the circulation of the individual vortices and Fp 

is the magnitude of the force exerted by the fluid on the bounding wall. 
The center of gravity of a vortex system is obtained from the relations 


y = 







z = 





( 5 ) 


for a y, z coordinate system. 
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Theorem 2: If a vortex system produces no moment with respect to 

its center of gravity the moment of inertia of the system will remain 
constant in value for all time, 

^ ^i r i^ = Constant (6) 


Here r^ is the distance of the vortex P^ from the system center 
of gravity. 

These theorems are written for a system of discrete vortices. For 
the case where the system is continuous (i.e, a vortex sheet ), the sum- 
mation V is replaced hy a line integration along the sheet. 


Deformation of the wake strip .- These results may now he applied 
to determine the deformation of the initially flat strip of wake taken 
from the vortex sheet immediately behind the trailing edge . A y, z coor- 
dinate system as shown in figure 7 is used in which the vortex strip 
initially lies along the Y-axis with the origin located at the center 
of the strip. Since the two halves of the strip are mirror images^ 
the plane y = 0 may be taken as a solid boundary and only the half 

0^ y ^ 

* 2 


need be considered. 


The intensity of the strip is initially 

v2 


ay 


1 - (fz 


2] -1 / 2 


and the center of gravity of the strip is 

ob/2 

^0 

y = - 


ar . 

y — dy 

<iy 


>b/2 


0 


dT y, 
— dy 

tfy 


z = 0 


itb 

8 


( 7 ) 
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The strip will initially tend to move downward as a flat sheet with a 

p 

velocity given by w = However, a flat vortex sheet is inherently 

unstable and begins to curl up at the edges as it travels downward. 
This instability may be physically explained in the following qualita- 
tive manner. The flow field about the flat sheet is identical to that 
about a translating plate where very low pressures occur at the edges 
due to the very high velocities. The fluid sheet is unable to resist 
these suction pressures at the edges and under the pressure gradient, 
rotational fluid from the interior flows into the tip region thereby 
increasing the circulation about the tip. (The vortex "sheet," it 
should be noted, is the idealized form for a finite thickness layer of 
rotational fluid. ) The initial vorticity distribution is thus altered 
and the velocity field produced is such that it causes the sheet to 
begin winding about the tip core. The rolling- up process during the 
very first moments of deformation have been treated in detail by Kaden 
(ref. 2). A numerical investigation of the rolling- up process using a 
row of discrete vortices to represent the vortex sheet has also been 
performed by Westwater (ref. 3). Both of these (two-dimensional) 


treatments indicate that the semisheet ^0 % y ^ coils up into a 
spiral of the form shown in figure 8. 


As the strip rolls up and moves downward 


moves along the line y = ^b, as predicted by 

8 

r_ 

velocity of the center of gravity is w = 


2b 


plete the velocity is w = 2 


?r 


K 


its center of gravity 
theorem 1. The initial 

and after rollup is com- 


Considering again the initially flat strip of figure 7j> the total 
circulation associated with a segment y^_ = y = ^ using equation (2), 
is 

^ ■ / 7 2 f * - r ° / m ¥) 2 (8> 

■ y i 

For convenience in subsequent derivations, y is now written in terms 
of the eccentric angle cp. 


y = | cos cp 


( 9 ) 
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Then, 


— = -r — cot cp 
dy ° b 

(10) 

r ( y l) = r o sin ^1 

(ID 


The center of gravity of the segment y 1 = y = is 


dT 


- F7^ Do y - * ■ 


r ( y l) ^b/2 dy **■ sin 'Pi 


K + \ sin ^i) 


z x = 0 


( 12 ) 


With respect to the origin y = 0, the moment of inertia of this seg- 
ment is 

J o( y l) = / y 2 ^ & = r o(|) 2sln 'Pl^ - j sin^cpp) ( 13 ) 

Transferring this moment of inertia to the center of gravity of the 
segment ^ y = ^ gives 


J ( y l) = J o( y l) - r ( y l) y l 2 


= r\ 


(I) 2sin M 1 - ? sln S - (t) 2 + 1 sln ^ 


(14) 


Now as the rollup proceeds, the strip coils up in a nearly circular form, 
as shown by the results of Kaden (ref. 2) and Westwater (ref. 3 ) j there- 
fore, the segment y^ = y = ^ coils into a spiral contained within a 

circle of radius r^. As the circular core grows, the circulation of 
the enclosed vorticity is given by 

r (n) - r (n) 


(15) 
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Similarly the moment of inertia of the segment after coiling is equal 
to its original value by theorem 2, 

J ( r i) = J (/i) ( l6) 

Thus when r increases by dr, y must decrease by dy (whence cp 
increases by dcp) and the corresponding increments in core circulation 
and moment of inertia become 


dT 

dr 

d T 

dr = — dcp = T 0 cos cp dcp 
dcp 

(17) 


dr 2.A dJ 

— r^dr = — dcp 

( 18 ) 


dr dcp 


Differentiating equation (l4) gives 


= r (1] cos^cp + - c . ° . s [cp + i sin 2cp 
dc P b sin 2 cp V 2 


2 sin cp 


^cp + i sin 2cpJ(l + cos 2cp) 


(19) 


Substitution of equations (17) and (19) into equation (l8) yields 


L 

1 

1 

5 

6 


4r£ 

b 2 


COS 


2cp 


4 sin^cp 


cp + ^ s i n 


ap) : 


cos cp 
sin cp 



(20) 


Since 


cp and y are related by the expression 


cos cp = - y. 


equa- 


tion (20) denotes the relation between y and r. A plot of — 

b 

against — as given by equation (20) is shown in figure 9- For 
b 

small values of cp (cp -» 0) equation (20) can be linearized by a 
limiting process to give 


2r 

b 


2 

3 


sin^cp 


( 21 ) 
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or in terms of y 



this becomes after a similar limiting process 


r 



( 22 ) 


A small segment of the strip near the tip is thus coiled within a circle 
of radius equal to 2/3 the length of the segment. Equation (22) has 
been plotted in figure 9 where it is seen that it closely approximates 
equation (20) through most of the range. Hence equation (22) may, with 
sufficient accuracy, be used to obtain the final core properties since 
over 95 percent of the total vorticity is contained within the core 
region where the agreement is quite good. Also shown in figure 9 is a 
p 2r 

plot of — against — . Figure 9 is taken from the results of 

r Q b 

reference 1. 


These results lead to the picture of two final vortex cores of 


radius r Q = — spaced a distance 


ix b 

4 


apart with the vorticity heavily 


concentrated near the centers, as shown in figure 10. 


Final wake form .- The deformation of the vortex strip is now traced 
as it passes downstream by using the axis system of figure 5 with the 
assumption that the deformation follows the relations just established. 
Everywhere outside the developing cores the x- component of velocity is 
taken to be V cos 5, while the y- and z- components are v and 
w - V sin 5, respectively. When the core has grown to a radius r 
the circulation of the enclosed vorticity is r(r). Since the cores 
are nearly circular in form, especially for small values of r, the 


tangential velocity v,p 


at the core surface 




+ (w - V sin 5) 


is assumed uniform and is therefore obtainable from the circulation. 


V T 


(r) 


r(r) 

2itr 


(23) 


The core grows by the addition of vortex filaments from the sheet to 
its outer surface. The filaments add along the streamlines which are 
determined by the velocity components V cos 6 and v,p to form con- 
centric vortex helices of intensity — . (See fig. 11.) This winding 

dr 

process continues until the original sheet has been consumed. 

By use of equation (22) relating the core properties to those of 
the initial sheet, the circulation, intensity, and tangential velocity 
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distributions within the fully developed cores can be determined as 
pure functions of the radius. The core circulation is 


r(r) = r D 


1 - (b - 5r) 2 
b 2 


1/2 


J 


(24) 


the intensity is 


^ = r ^ (t - 5r) 

dr 0 b 2 


1 - -i- (b - 3r Y 


- 1/2 


(25) 


and, by using equation (23), the tangential velocity is 


V T " 2itr 


1 - (b - 3r) £ 
b 2 


1/2 


(26) 


It is evident from equations (25) and (26) that both — and v T 

dr 

become infinite at r = 0. These are also properties of a mathematical 
point vortex. Unlike the point vortex, however, the velocity distribu- 
tion given by equation (26) leads to a finite kinetic energy for the 
core flow, as will be shown later. 


Axial velocity field u(r) .- The preceding development led to the 
picture of the cores as a superposition of concentric vortex helices, 
each helix being a vortex filament from the original sheet. A helix of 


radius 


r 


has an intensity 


dT 

dr 


and induces a constant longitudinal 


velocity du across its interior. As shown in appendix A the differ- 
ential induced velocity due to this elemental helix is 


du = n(r) — dr (27) 

dr 

where n(r) is the number of vortex turns per unit length of the core 
and 


L 

1 

1 

5 

6 


v T (r) 

2nr V cos 5 + u(r) 
_ — 


n(r) = 


(28) 



17 


At the surface of the developing core the longitudinal velocity is 
V cos 6, since u = 0 ; therefore, n(r) is given by 

v T (r) 

n(r) = 

2 itrV cos 5 

However, as more vortex lines are added, the longitudinal velocity u(r) 
within the core tends to change, thus changing n(r) of the filaments 
within the core. After the cores are fully developed, n(r) is given 
by equation (28). The value of u at a given radius is obtained by 
summing the velocities induced by the helices exterior to that radius, 


ob/5 

u(r) = / n(rj) •— dr) 

J v <17 


dT 


(29) 


where t) denotes a dummy variable of integration. For the region 
— = r = 0, using equations (26) and (28) gives 

5 


n(r) 


iv 


(2nr) 2 [v 


cos 5 + u 


(r)| 


1 - 1 (h - 3r) 2 ! / 
b 2 i 


(30) 


therefore. 


“(=-) - -f-5 r 0 2 

4 * 2 b 2 


'b /3 


r |v cos 5 + u(r)) ^2 7 j 


b 


(3D 


This integral equation has the equivalent differential form 


2 

(V cos 5 + u) ^ - (± - f| 

^ 4 A 2 \r 2 


The solution of this equation, subject to the condition uj- 2 J = gives 

the axial velocity field. Although the differential, equation is non- 
linear, it is readily integrated after separating the variables. Thus, 



(V cos 5 + u’ )du' 
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Expansion of this equation yields 


1 

2 


u 


+ uV cos 


5 


9 r p 2 


lo S e r 



1 


and. since 



cos 


2 5) 


the axial velocity field becomes, as a function of core radius and wake 
inclination angle. 


u 



V 2 cos 2 5 + 




V cos 6 (52) 


or 


where 



l)fl(r) - 1 


fl(r) 



1 - 



(35) 


This relation (eq. ( 52 )) predicts an infinite value for u at 
r = 0 , but unlike the velocity field v T , the kinetic energy of the 

u(r) field is infinite. In order to obtain a more useful picture of 
the axial velocity field, divide both sides of equation ( 33 ) by V cos 6 
to obtain 


u 

V cos S 



(tan 2 s)f 1 (r) - 


1 


( 5*0 


which gives the ratio of the axial velocity u to the wake- direction 
component of the free- stream velocity. A plot of f^(r) against r/b 

is presented in figure 12. Now for small values of 5, say 5 = 5°> 

tan 2 S = 0.008. For most of the core area, say r = 0.1b (this includes 
91 percent of the total core area), f^r) is small, reaching the value 

of 1.1 at r = 0.1b. Thus for small wake angles the axial velocity u 
is of negligible magnitude over most of the core area. Very near the 
center, however, the ratio attains very large values, becoming infinite 
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at the core centers, even for values of 5 approaching 0. It is 
obvious from physical considerations that the large core velocities 
cannot exist, and the reason for their prediction by equation (34) can 
be directly traced to the assumptions made in the derivation of this 
result. 

In the derivation of equation (29) (see appendix A), it was assumed 
that n was sufficiently great that each vortex turn could be considered 
to lie in a plane. This assumption is equivalent to a pitch approaching 
zero for all helices. However, an analysis of equation (30) with u = 0 
shows that even for appreciable values of 6 (5 ^ 30°), n is small 

(n « l) for the outer helices (r near b/3). Hence, the planar con- 
dition is not approached and due to the inclination of the induced- 
velocity vector with respect to the core axis the potential magnitude 
of the axial velocity over the core interior will not be realized. 
Additionally, these same effects on interior helices may produce nega- 
tive axial velocities which can counteract the positive induced veloc- 
ities of the outer helices. 


In the derivation of equation (26) the rollup of an infinitely 
thin vortex sheet was considered, so that the vorticity at the tips 



was infinite. 


In reality, the vortex- sheet wake consists of 


a layer of rotational fluid which attains appreciable thickness near 
the tip (ref. 4 and fig. 13). Since the vorticity of the layer is 
spread over a finite volume of fluid the vorticity density at the tip 
cannot be infinite, for a finite circulation in a circuit surrounding 
the tip region. Surveys show that near the center of the rolled- up 
vortex cores the fluid possesses a large but constant vorticity and 
the fluid there must rotate as a solid. Within this small subcore, 

Vrp — > 0 as r 0, so that n again becomes increasingly smaller as 


the center is approached. Thus, no infinite axial velocities can 
actually occur. 


These effects associated with a more accurate representation of 
the physical flow regime as it actually exists, offer a sufficient 
explanation as to the invalidity of equation (34) very near the core 
centers . 


It is of course possible to extend the analysis to a more refined 
consideration of the physical flow regime so as to include the effects 
of the subcore on the axial velocity field. The primary problem then 
is to establish the size of the subcores. However, due to the com- 
plexity of such a treatment, and in view of the fact that u is very 
small, even for relatively large values of 5 (5 = 3®°)> as is con- 
firmed by experimental evidence (ref. 4), such an analysis is not 
justified for the purposes of the present investigation. Therefore, 
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it is now assumed, to the degree of accuracy of the foregoing assump- 
tions, that the axial velocity field is sufficiently small so that its 
effects may be neglected in the calculation of the airfoil forces. 


Rate of Deformation of the Vortex Sheet 

The preceding results concerning the deformation of the vortex 
sheet can be used to determine the degree to which the sheet has rolled 
up at any given downstream station ("downstream" referring to distance 
along the X-axis in fig. 5)* Considering a cross section through a 
developing core, such as shown in figure 14 where v T is the tangential 
velocity at the core surface (assumed uniform around the periphery), the 
rate at which the sheet adds to the core is given by 


= 

dt 


V T 


(35) 


where dy is an elemental length of the sheet. The relation between y 
and r as previously developed is (eq. (22)) 


y 



whence 



(36) 


The time required for the core to develop to a radius r is therefore 


t- / i-^4 

J 0 2 v t(W 

Substitution of v T from equation (26) yields 


t = 


nVb sin 5 


o L 


1 - 4 (b - 3 T)) 2 


- 1/2 


dr] 


(37) 


The downstream distance x corresponding to this time is 


x = Vt cos S = — cot 5 
rtb 


1 - 4 (b - 3 t 


- 1/2 


dn (38) 
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The rolling-up process is complete when r = -4, and the time and 

5 

distance for complete development are obtained by setting the upper 
limit equal to b/3 in the common integral of equations (37) and (38). 

Using the general result 

r ~2 


(a 2 - p 2 ) 


1/2 


dp 




.-1 P 


a - p + — sin"- 1 - — + Constant 
2 a 


and associating p with tj (p = rj 1 / 2 ), the common integral of equa- 
tions (37) and (38) can be evaluated to yield the results 


t = 


4b 


nV sin 5 


- - \ - \ - - - + - sin -1 1/ - 
2 V b V 3 b 3 V 2 b 


(59) 


4b 4. c 

x = — cot 8 
n 


- 1 iff (/FI + \ Bin ' 1 /li 


(40) 


for any value of r. 

For complete rollup of the elliptical sheet, r 


— , therefore, 


r = 0.1211 ^ cos 5 


(4l) 


d = 0.1211b cot 5 


(42) 


For small angles. 


Since for an elliptical wing 


& = JL £2 

*2 Vb 


6 = — — 

n 5 A 


equations (4l) and (42) may be written in the following form: 


x =0.93 | A 


(45) 
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d = 0.93b jjL (44) 

This form is useful for comparison with Kaden's results (ref. 2). 


DETERMINATION OF THE AERODYNAMIC FORCES 


The lifting airfoil experiences a resultant aerodynamic force F 
which is composed of the lift and induced drag components taken nor mal 
and parallel, respectively, to the free- stream velocity. In the present 

system of axes (fig. 5) the aerodynamic force F is resolved into com- 
ponents F z and F x as shown in figure 15 where the positive force 

directions are indicated by the vectors. The transformation equations 
for the two force systems are, from figure 15, 

L = F z cos 6 - F x sin 5 (45) 

Di = F z sin 6 + F x cos 6 (46) 

Determination of any two force components allows the calculation of the 
resultant aerodynamic force F; therefore, F z and D-^ will be deter- 
mined for this purpose. 

Determination of F z 

Far behind the airfoil the wake consists of two vortex cores con- 
taining the coiled vortex sheet. Within the vortex cores the motion 
is irrotational except for the surface of velocity discontinuity repre- 
senting the coiled sheet. This wake has been created from air which 
was initially at rest (considering the flight of the wing through still 

air). The creation of the motion associated with a unit length of the 

— > 

wake has required the application of a definite vector impulse -I 1 by 
the airfoil. From the viewpoint of a stationary airfoil in a free- 

stream flow of velocity V (steady flow), the instantaneous time rate 
of change of the total impulse of the wake is constant and equal to a 
component of the aerodynamic force acting on the control volume, 

(47) 
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may be resolved into its 


m 


^Fy' =0 because of the plane of symmetry of the wake j . 

The coiled vortex sheet may be considered, as previously mentioned, 
as being formed from a superposition of concentric helical vortex fila- 
ments of strength — dr. Corresponding filaments (one from each core) 

dr 

form a vortex pair which may be considered part of the boundary of an 
infinite vortex ring of which the bounding filament is helical in form. 

The total impulse of a unit length of the wake may therefore be obtained 
by summing the impulse contributions of all the vortex filaments com- 
prising the cores. 

For any vortex ring formed by a closed vortex filament of strength T 
the impulse acting on the boundary generating the motion from rest is 


Since the impulse is a vector quantity, il 
components, whence 


F ’ = 


dl; 

dt 




dt 


V - 0 


?= ff S P(?a- *b)^ 


(^ 9 ) 


where p0 is the impulsive pressure acting at a point on a surface S 
having the ring as its boundary, and 0 a and 0^ refer to the potential 

values at neighboring points on either side of the surface (fig. 16). 

On passing through the surface the potential must suffer a jump equal 
to T whence 0 a - 0^ = F, and 


I = IT pfn dS 


( 50 ) 


where n dS = dS, gives the resultant linear impulse to generate the 

— > 

ring. The value of I is independent of the form of the surface S 
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taken, by applying Stokes' theorem. Since P is constant by Kelvin's 

q’.d? = 0 , the force required to increase 

the impulse of the ring is 

= p r JL ff n dS (51) 

dt dt JJ 

The impulse of a helical vortex pair located at radius r within 
the core is now determined. The axes of the helices are spaced kb/4 
apart as shown in figure 17- Corresponding points on the two filaments 
may be connected by a straight line parallel to the Y-axis so that a 
surface may be constructed, all of whose generators are normal to the 
X- and Zr-axes and connect corresponding points on the two helices as 
shown in figure l8. If this area is denoted by S, the impulse per 
unit length of the pair is 


circulation theorem 


m 


dl ' 



p — drn dS 
dr 


(52) 


The z- component of this impulse is 


dl 


t 

z 




— > — 

n.z dS 


(53) 


— ^ ^ 

where z is a unit vector parallel to the Z~axis. However n.z dS is 

merely the projection of the area n dS on the X-Y plane (fig. 19 ), so 
that equation (53) becomes, considering the symmetry of the projection 
of a helix with regard to its axis. 


4 


dl z ' = p * b dr 
z > dr 


( 5*0 


and this gives the z- impulse per unit length of a helical- filament pair. 

From the results of the section "Axial velocity field u(r)," the 
longitudinal velocity u is negligible so that integration of equa- 
tion (54) yields the simple result 


_ nb ar 

■ F z = J„ P - V cos S - dr - p 


nb 




cos 6 


(55) 
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for the z-component of the aerodynamic force. This is exactly the same 
result as would he obtained by considering the cores as line vortices 
of strength T 0 . 

From figure l8 it is evident that there should exist also a force 
component F x ! , since from the geometry of the wavy impulse surface it 

can be seen that 



x.n dS < 0 


where x is a unit vector in the positive x- direction; that is, there 
is a net component of impulse surface area in the negative x- direction. 
This force F x f is clearly the force associated with the increase in 
longitudinal momentum due to u. Under the assumption that u is 
negligible, F x * =0. 

A second method for determining F z through a consideration of 

the hydrodynamic mass of the wake body is given in appendix B. The 
results, as may be expected, are again identical tc those for a vortex 
pair. 


Determination of D± 

The component Dj_ of the aerodynamic force is determined by an 

integration of the wake kinetic energy. The basic relation for this 
calculation is the power equation, 


^ = VDi = pV cos 
dt 



v2 + w2 
2 


dS + pV cos 5 



+ v 2 + 


dS 


+ P 



2 


u dS 


(56) 


where °°, 0 denotes an integration over a plane taken perpendicular to 
the cores and far downstream of the wing, of the area external to the 
vortex cores, and S c denotes integration over the core areas. If u 
is considered negligible, equation ( 56 ) becomes simply. 
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= p cos 5 


II 


(a) 

“'° v2 + w 2 


0>) 


dS + p cos 5 


S c v 2 , w 2 . . 

— - dS (57) 


These integrals are designated (a) and (h) as shown. 

Integral (a) is simply evaluated by considering the two-dimensional 
velocity field due to the vortex pair. The complex potential is (fig. 20) 


W = iK log e - iK log e R2 


(58) 


r o R 1 

where k = — . The stream function of the motion is i - k log so 

2 it c R2 

that 


cos 

2 


b If ’ (t 2 + ^ = f oos If ’ ( w l v !■) 


v dy dz 


P 

= — cos 
2 


b If ’ fe (m,,) - 1 


dy dz 

( 59 ) 


since the flow is everywhere irrotational. By Stokes' theorem and 
equation (59) 


p 

| cos 


= | cos 8 /%( v dy + w dz) 


nC 


( 60 ) 


where the contour C denotes the core boundaries and the connecting 
line as shown in figure 20. The integration is carried out in the 
positive direction with the enclosed area to the left. Since 
(vdy+vdz) = vrj\ ds, where vrp is the tangential speed at the core 

boundary, equation (60) becomes 


cos 


HP °°,0 v p C^_ 

6 JJ (v^ + w^jdS = -p cos 8 (jf) KVfj, log e U- ds (6l) 


ht H H tf\VD 
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■where denotes counterclockwise integration around the right core 

only. On this boundary Ri = ^ so that 

5 


■§ 


Cl 


kv t log e r x ds = - — log e - 


( 62 ) 


r l b 

rgb’’ 


In the expression log e ~ = log e b is cancelled before inte- 


( 

gration . J For we have 

p Ci o 2rr p 

K(log e r 2 )v T ds = - J V T 2rt l0ge r 2^ 9 


)de (63) 


From the geometrical relationships of figure 21, the following functions 
of 0 are obtained: 


r 


v t ( 6) p- ph) 


cos 


. • -1 / Sin 0 

0 - sin x - — 7 — r- 

l 3r 2 (0) 


r o(e) = 


^ + i COS 0 1 + 

4 3 




1/2 


Graphical integration of equation (63) then yields 


(64) 


r c i 


(j) K ( lo g e r 2) v T ds = - 6.432 


12jt £ 


(65) 


and integral (a) of equation (57) becomes 


00 J 0 X 

~ cos 6 /j (v 2 + w 2 )dS = 0.1205P cos 6 r o 2 


(66) 


Integral (b) of equation (57) is to be evaluated over both vortex 
cores. If the geometrical relations of figure 22 are used, the core 
velocity components are 
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v = 


2* 


r sin 9 


* 2 p 

— b + r cos 0 'i + (r sin 9)^ 

4 / 


- ^ (b - 5r)2 


w = 


2it 


— b + r cos 9 


* '. 2 

— b + r cos 0i + (r sin 0) 

4 J 


cos 9 
r 


1 - 4 (b - 5 r )‘ 

D 


l (67) 


The velocity v T within the cores, as given by equation (26), 

becomes infinite at r = 0. The kinetic energy of the core however 
approaches a finite value. If a very small region near the center of 
the core where viji becomes large is considered, the kinetic energy of 

this region can be obtained by summing the rotational and the transla- 
tional energies. Thus, 


dE 


' = P v T^ 2ltr ^ + V 2 sin 2 5 i p2n:r dr 


where 


p 2 , 

v 2 = 12_ _L 

T 4«2 r 2 


1 - 4 (* " 3r) : 


and 


E’ = < 


P P r 


4 it 


I n 

0 71 


4 (b - 3T Y 
b 2 


dr) + 4 4 P r < 

b 2 


_1 

4it 


U £ - 2 £) + 2. £ 

lb 2 b 2/ b 2 


P r o" 


( 68 ) 


L 

1 

1 

5 

6 


and the kinetic energy per unit length of the core is finite, even 
though the rotational velocity is infinite at r = 0. 

Now arbitrarily choosing a small value for r, say r = 0.001b, 
integral (b) can be evaluated by using equation (68) to obtain the 
induced drag corresponding to the energy in the region 0 % r = 0.001b 
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= 9.1098 sin 2 6 cos 5pV 2 b 2 


(69) 


Aerodynamic Forces 

Solving equation (46) for F x and substitution into equation (45) 
gives the lift force, 

L = F z sec 5 - Dj, tan 5 ( 70 ) 

From equation (55)> after substituting for P , 

F z = 3-8757 sin S cos SpV 2 b 2 ( 71 ) 

From equations ( 69 ), ( 70 ), and ( 71 ), the final expressions for the lift 
and drag forces are 

L = pV 2 b 2 (3.8757 sin 5 - 9 .IO 98 sin ^5 ) 

= pV 2 b 2 (9.1098 sin 2 5 cos 5) 


L 

1 

1 

5 

6 


Since — = A, the corresponding coefficients may be expressed as 


C L = 2A(3-8757 sin S - 9.1098 sin ^6 ) 
= 2A(9.1098 sin 2 5 cos 5) 


(72) 


It follows from equations (72) that the forces are directly proportional 
to the aspect ratio A. 

The results of equations (72), obtained by using 5 as a parameter, 

c c 

are plotted in figure 23 in the form — against — The linearized 

A A 


Cl 

lifting-line result — = \ln 


c Di 


is included for comparison. 
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and by using the 
to the energy in 


following integral for the induced drag corresponding 
the region 0.001b = r = 

3 


p cos 5 


2 rb/3 p2jt 
4rt 2 ^O.OOlb^O 


r 


2 


sin 0 


1 - ^ (b - 3r)‘ 
b^ 


r sin 0 

^ b 4- r cos 0^ + (r sin 0)' 


b 4- r cos 0 


- £2^1 J 1 - i (b - 3r) 2 


b 4- r cos 0^) 2 4- (r sin 0)* 
ir dr d0 


J 


This integral may be evaluated numerically for each core yielding the 
total value (sum of both cores) 


5 . 0056p cos 5 — 

2* 2 

The integrations were performed by digital computer using Simpson* s 
rule. Evaluation of equation (68) for the region 0 = r = 0.001b 
over both cores yields for the induced drag 

- (0.0060 - 0.0000045) + — Q0 - 00Q -- 

3 71:2 

This subcore value is an extremely small percentage of the total kinetic 
energy of the cores and an even smaller percentage of the total flow 
kinetic energy. Hence it is immaterial whether the velocity field of 
equation (26) or that due to a subcore of constant vorticity (as is 
the actual case) is used to determine the kinetic energy of the flow 
near the core center, since the value is in either case negligible. 

This is true, of course, only if the radius of the constant vorticity 
subcore is small, and experiment shows this to be the case (ref. 4). 

Substitution of the final values for integrals (a) and (b) into 

2 

equation (57) an^ replacing T Q with — Vb sin 5 gives the induced 
drag as a function of 5 as follows: 


p cos 5P r 


jt 
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DISCUSSION OF RESULTS 


One of the most important results to be noted from the normalized 
drag polar of figure 23 is that a maximum lift coefficient exists and 
this maximum coefficient is proportional to the wing aspect ratio. From 
equations ( 72 ), the maximum lift coefficient is 

c L,max = 1 *9^A (73) 

and this occurs for a wake angle 5 = 22.1°. Thus C^max obtained 
when r Q = 1.86Vb for an elliptical circulation loading. 

The relation for Ct^ may (eq. (73)) requires some explanation 
since it predicts unrealistic values of lift coefficient for wings of 
large aspect ratio. For example, for A = 15, CL,max = 29*1. It is 

clear that for wings of large aspect ratio. Cl, max is determined by 

some effect other than the wake induced velocity. For normal wings 
such as considered herein, Cp is set by the physical ability of 

the section to produce circulation. If a Joukowski profile, as obtained 
by conformal transformation of a circle of radius a is considered, the 
theory shows that for the section, 

C L = 8rt | sin 

where c is the chord of the profile and ct^ is the absolute angle of 
attack of the zero lift line. This angle a a determines the circula- 
tion around the section, since the Kutta condition requires that the 
point on the circle corresponding to the profile trailing edge be a 
stagnation point. The maximum value of a/c is l/2, so that for this 
case 


- k-n sin a a 


and 


< -'L,max “ 

for a & = 7L, When <x & = the front and rear stagnation points for 

both the circle and the profile have coincided. For the profile, both 
stagnation points then lie on the trailing edge, and any further pro- 
duction of circulation is impossible from a physical standpoint. If 
all profiles are assumed to behave in a manner similar to the Joukowski 
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profile, the limit C L ^ max = 4 jt = 12.6 may be set for all profiles. 

Nov assuming that a maximum section lift coefficient of 12.6 can be 
attained by each section of the three-dimensional wing, the wing lift 
coefficient cannot exceed 12.6. Hence, regardless of aspect ratio, 
the absolute Cl, max the cannot exceed 12.6. This leads to 

the conclusion that only wings having A = 6.5 will have Cj^ max 


limited by wake effects. Here 



where 


&w - 



These results are summarized in figure 24 which is a plot of maxi- 
mum circulation lift coefficient against aspect ratio. Below an aspect 
of 6.5 CLjniax gi ven b y the relation 


^L,max " 

Above this value of A, max is equal to 12.6. This curve may be 

considered the boundary of all possible circulation lift coefficients. 
(See fig. 24.) 

A brief analysis of the physical meaning of Cp/A is of interest. 

In the derivations only an elliptical distribution of circulation is 
considered to exist for each value of r o , without regard to the physi- 
cal means for obtaining such a distribution. In practice any number of 
wing configurations producing a given lift coefficient and having 
elliptical circulation can be had merely by a suitable combination of 
chord, camber, and twist distributions, for the same span. Hence the 

b 2 

aspect ratio A = — can scarcely have any unique physical meaning, 

^w 

since it depends upon the particular chosen for its definition. 

The value of S w to be used in computing the value of A for use with 

figure 23 is therefore arbitrary, so long as the same value is used to 
define the lift coefficient for the wing. This follows from the fact 
that the normalized lift coefficient is independent of wing area: 


1 = _L_ fh[ = _1_ 

A <l S w b 2 qb 2 


L' 

1 

1 

5 

6 


I 
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1 
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Thus, figure 23 may he used as the drag polar for wings which appreci- 
ably alter their geometry and area for various regions of the lift 
range, provided of course that the circulation distribution at each 
lift condition approximates the elliptical form. The problem of 
designing a wing which will possess an elliptical or any other speci- 
fied distribution of circulation at a given necessitates a knowl- 

edge of the downwash distribution at the wing and, as previously men- 
tioned, this information cannot be furnished by momentum considerations. 


is shown in fig- 


= 1.1. Thus use of 


c L ( CpiA 1 / 2 

The Prandtl lifting- line result — = h jt — — \ 
ure 23 for comparison with equations (72). No appreciable deviation of 

A 
*A 

appears valid up to the values of lift coefficient equal to the aspect 
ratio of the wing. This result is also summarized in figure 2k which 
shows a plot of the equation = 1.1A. For combinations of Cp 


the linear and nonlinear relations occurs below 

.2 

the linear- theory result for calculation of the induced drag 


and A below this line, lifting-line theory is valid. Above this line 
the nonlinear relation of equations (72) is applicable. 


Unfortunately, experimental data at sufficiently high values of 
Cp/A for comparison with the nonlinear region of the theoretical 

results do not exist. However considerable data are available for 
boundary- layer- control wings of various aspect ratios and planforms 
operating at moderate values of C-^, and experimental results (refs. 5 

to 15) for several wings are compared with the theory in figure 25. 

As shown in the figure, the experimental points represent a range of 
aspect ratios, planforms, and boundary- layer- control systems. In all 
cases the maximum lift coefficient of the wing is limited by stall, 
and the points shown represent values just before the stall of the 
wing. The data correlate well with the normalized polar, even for the 
rectangular planform and wing-body configurations. The data for the 
wing of aspect ratio 20 correlate especially well, as might be expected. 
Reasonable agreement occurs even for wings with appreciable amounts of 
sweepback. This agreement may be considered as experimental verifica- 
tion of the theoretical prediction that the lifting- line relations are 
valid even for high lift coefficients provided the aspect ratio is suf- 
ficiently large % 1.1A). 


C 


L 


An additional point of interest in figure 23 is the relation 


= C 


Di’ 




= 1.90, 


and if the total drag of the wing 
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is assumed to equal the induced drag the thrust required for flight is 
equal to the lift. Thus, an aircraft capable of operating at C L = I.9OA 

must have sufficient thrust to take off vertically. This result may be 
used as a rough demarcation line between STOL and VTOL capabilities of 
a given STOL aircraft. 

The drag polar of figure 25 as established by this investigation 
is similar to that arrived at by Helmbold (ref. 14) using a different 
procedure. The agreement is due, however, to the fact that the axial 
velocity field is assumed negligible in the present treatment. In the 
derivations of reference 14, no restrictions were placed on the induced 
axial velocity so that the agreement with the present results indicates 
that the assumptions of reference l4 reduce it to a two-dimensional 
treatment of the problem. 

The result given in equation (42) predicting downstream distance 
at which the vortex sheet is fully rolled up may be compared with the 
result of Kaden (ref. 2). Kaden's relation for an elliptical sheet, 

d = 0.28b — 

CL 

may be put in the form 

- = 0.04 cot & 
b 

Cr r 

for small values of 5, since — = — — . From equation (42), 

' A 2 Vb 

£ = 0.12 cot & 
b 

The present theory predicts a distance three times as large as does 
Kaden's result. However, Kaden's derivation is based upon the rollup 
of a semi- infinite vortex sheet, so that the question as to when the 
rollup is complete is indeterminate, making the coefficient 0.28 some- 
what arbitrary in nature. The result of the present theory is based 
on complete rollup of a finite semisheet. The apparent difference in 
the two results is practically negligible, however, as may be seen by a 
consideration of figure 9* Here it is seen that the greater part of 
the vorticity is concentrated near the center of the core so that, for 
example, when the core has grown to only 64 percent of its final radius 

= 0.44^ it already contains 92 percent of the vorticity. Using 
this radius, r = 0.22b (where equations (20) and (22) diverge 
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appreciably), as the basis for essentially complete rollup, the result 
from equation (40) is 

- = 0.05 cot 6 


which may be taken as the "essentially" complete rollup distance. This 
choice of radius is, however, purely arbitrary and the complete rollup 
distance is given by 


£ = 0.12 cot & 


That the complete rollup requires considerably more distance than 
the essentially complete rollup is shown in the experimental data of 
Muttray (fig. l8 of ref. 15 ) where the vortex cores move inward lat- 


erally toward the equilibrium line y 


-f b 


only very slowly. 


A knowledge of the rollup distance is primarily of interest in 
estimating the downwash field for longitudinal stability analyses. In 
lieu of accurate experimental data, a first-order approximation of the 
downwash field behind a highly loaded wing might be estimated from the 

foregoing results. For this purpose the relation — = 0.05 cot 5 is 

b 

plotted in figure 26 as a function of Cl/A by using equations (68). 
This plot shows for example, that the vortex sheet of a wing with 


, / C T \ 

aspect ratio 6 pulling a lift coefficient of 5 y— = 0.5 j will be 

essentially rolled up at a distance d = 0 . 76 b along the wake axis. 

If the tail location is such that it lies beyond d for. a given Cl /A, 


the velocity field may be estimated by use of equations ( 58 ) and ( 67 ). 
If the tail is located forward of d the velocity field might be 
roughly estimated by use of Kaden's results (ref. 2) for the lateral 
location of the cores, combined with equation (40) which determines 
the size of the cores. It is emphasized that this method is quite 
approximate, but it furnishes a means for first-order analysis of sta- 
bility effects at high lift coefficients. 


CONCLUDING REMARKS 


The objective of this paper has been the development of the com- 
plete drag polar for wings having an elliptical distribution of circu- 
lation and the objective has been sought through a detailed considera- 
tion of the wake produced. In the course of the derivations the 
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assumption was made that the axial velocity within the cores was negli- 
gible in its effect on the wing forces and an attempt was made to 
justify this from both theoretical and experimental considerations. In 
view of the other simplifying assumptions made in the derivations a 
more detailed analysis of the axial velocity field does not appear 
justified. 


The resulting theory predicts that the wake induced velocity limits 
the max imum lift coefficient to a value equal to times the wing 

aspect ratio A for wings with A ^ 6.5 and that airfoil section char- 
acteristics limit the lift coefficient to 12.6 for wings with A > 6.5. 
The theory also indicates that up to lift- coefficient values of 1.1A 
the results of the linear lifting- line theory are valid. Although 
experimental data for values of lift coefficient Cp sufficiently high 


to compare with the theory near 


/£l\ 

\ A/ 


are not available, existing 


max 


data for many current wing planforms using boundary- layer control cor- 
relate well with the theoretical results in the lower range of C l /A. 


Experimental data indicate that with the relatively small amounts of 
boundary- layer control currently being used on operational aircraft, 
the stall lift coefficient is still sufficiently low that linear theory 
is applicable. 


Finally, experimental evidence indicates that the theory predicts 
with reasonable accuracy the lift- drag relationship near stall Ci^ max 

of wings with widely varying planforms, flap deflections, degrees of 
sweep, and methods of boundary- layer control. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., November 14-, i960. 
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APPENDIX A 


THE AXIAL VELOCITY FIELD OF A VORTEX SOLENOID 


The vortex lines comprising the vortex sheet which emanates from 
the trailing edge of a lifting airfoil coil about the tip filament to 
form concentric helices as they pass downstream. By virtue of their 
helical form they act as a solenoid, inducing an axial velocity u 
within the vortex core. The following is a brief derivation of the 
velocity profile of such a vortex solenoid. 


Magnitude of u Due to a Single Vortex Ring of Radius r 

Far downstream the solenoid can be considered as having infinite 
length and being circular in form. The number of turns per unit length 
of each helix is assumed sufficient that each turn can be considered to 
lie in a plane normal to the axis. From induction theory, the axial 
velocity induced at its center by a ring of radius r is 


^ u center 


d r 
4jrr 2 



d LT 
2r 


(Al) 


where dT is the strength of the vortex filament and 9 is the angle 
as shown in figure 27. For any point P on the axis of the helix, 
located a distance x from the ring, the axial velocity induced by the 
ring is 


Au- 


P, axis 


dTr 


, 2 * 


4jt(x 2 + r 2 


\3/2 


r d.9 = 


r 2 ar 


J 2 7 2 ^ 3/2 

2\x + r ) 


(A2) 


Magnitude of u Due to a Single Infinite Helix of Radius r 

If n is the number of turns or vortex rings per unit length of 
the helix of strength dF and radius r, then the circulation per unit 
length is n dT. The circulation about a ring of width dx is n dr dx 
(fig. 27). From equation (A2) this ring will induce at a point P on 
the axis the velocity 


^ U P, axis 


P P 

nr dP dx _ nr dT dx 

2(x 2 + r^f^ 2 r3 
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where R is defined in sketch at "bottom of figure 27* For the velocity 
at P due to a length l of the helix (using bottom sketch of fig. 27), 


Au- 


P, axis * ™ ® 9 1 - =° s 0 2) 


l0 2 


n _di\ 

2 


since 


n dT n dr R d9 „„.„2 0 n dr n 

r dx = sin 9 = sin 9 d9 


2R3 


2R sin 9 


When the helix is infinite, ®l = ®2 = ^ and 


Aup = n dF 


(A3) 


Velocity Profile u(r) for Entire Vortex Core 


A helix 
strength dr 


of radius r having n(r) turns per unit length and 

= — dr contributes an element of axial velocity of 
dr 


du = n(r)4^ dr 
dr 

at all points on the axis of the helix. Since the helix is infinitely 
long this velocity is constant over the interior. By summing all the 
helices comprising the core, the axial velocity profile is 


u(r) 


k 

lim ) Au^ 
k -» 00 i =1 



n(T))~ dT l 

dr] 


(Alt) 


where T) is a dummy variable of integration, r Q is the outer radius 
of the vortex core, and k is an integer. 
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APPENDIX B 


ALTERNATE METHOD FOR CALCULATING F z 


The stable wake far downstream from the wing consists of two vortex 


cores of diameter 


2b 


with centers 


-h 


apart moving normal to them- 


selves with a velocity V sin & (fig. 10). The relative closeness of 
the cores results in the formation of two distinct systems of streamlines 
relative to the cores, as will he shown. One system is closed - that 
is, the streamlines form closed curves about the cores - whereas, the 
other consists of open streamlines - that is, beginning and ending at 
infinity (fig. 28). The result is that the vortex pair entraps the 
fluid within the closed system and carries it along in its translation 
with a velocity V sin 5. The force exerted by the airfoil to give this 
mass of air the translational velocity V sin 5 is clearly the aero- 
dynamic force F z . This force consists of two parts: that associated 

with the translational momentum of the closed body of fluid and that 
necessary to accelerate the body of fluid to the velocity V sin &. 
Determination of the time rate of change of the z- impulse associated 
with the formation of this wake will, as before, give the force F z . 


The calculation is begun by first determining the boundary of the 
closed body, that is, the enclosing streamline, from which the sec- 
tional area may be established. Since each core has a strength k q 
the complex potential of the steady flow is, in terms of the complex 
variable z 1 = y + iz, 




l°g e (z ' 



i 

+ (V sin 8)e 



(Bl) 


whence 


* = £ log * 


[(y - y x ) g * s 2 ] 1/a 

[(y - y 2 f ♦ 4' 2 


- Vy sin 5 


(B2) 


Here y and z are Cartesian coordinates. The bounding streamline is 

2 

given by t = 0, and since r o = ^- Vb sin 8, the equation of the boundary 
becomes, after some transformations, 


l+o 



This implicit equation can be solved graphically to obtain the boundary 
streamline. The section of the wake is an oval with minor axis of 

0.926b and major axis of 1.659b* with an area of 1.760b 2 . As is evident 
from equation (B3), the cross-sectional area is independent of both P 0 

and core radius. 

The wake may thus be considered as a fluid cylinder of oval section 
moving normal to its axis with velocity of V sin 5. This cylinder 
carries a momentum per unit length of 

pA'V sin 5 (B4) 

where A' = 1 . 760 b 2 , the sectional area of the cylinder. The part of 
the force F z due to the cylinder momentum is obtained from the time 
rate of change of the wake z-momentum due to the creation of the cylinder, 


(F z )^ = pV 2 sin 8 cos 8A' = 1. 760pV^b^sin 8 cos 8 (B5) 

The second part of the force F z is that associated with the flow 
momentum external to the cylinder. The fluid motion external to the 
fluid cylinder is identical to that associated with the motion of a 
solid cylinder of the same shape. As predicted by hydrodynamic theory, 
the application of a definite impulse is required to accelerate an 
immersed body from rest to a steady velocity U in a fluid also at rest 
initially, and this impulse is greater than that needed for equal accel- 
eration in a vacuum. Thus the presence of the fluid increases the 
effective mass of the body by an amount M', the hydrodynamic mass, 
which is equal to the drift mass of the particular body. 

The power expended in generating a unit length of the exterior 
flow about a solid cylinder of any sectional shape moving normal to its 
axis is 


— = F'U = M'U — (B6) 

dt dt 

where F' is the part of the total force which must be applied to the 
body to overcome its acceleration drag. Converting equation (b 6) to 
an impulse gives 
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F'dt = M'U = I' 


4 l 


(B 7 ) 

where I' is the impulse per unit length of cylinder to generate the 
external motion from rest. The determination of the second part of F z 
now reduces to a calculation of M' for the oval cylinder of the wake. 

The oval is assumed to "be a true ellipse, since it varies only 
slightly from an ellipse with similar minor and major axes. The kinetic 
energy in the flow exterior to an elliptical cylinder with major axis 2a 
and minor axis 2 b, moving v/ith the steady velocity U normal to its 
semiaxis a is, per unit length, 

E = | pnl^a 2 (B8) 

Using equation (b 6) yields 

M’ = pita 2 (B 9 ) 

for the cylinder drift mass. In the present case a = 0 . 819b and 
U = V sin 8, and I' is obtained from equation (B 7 ) as 


I' = ( 0 . 8 l 9 ) 2 p*b 2 V sin 5 


Then 


(F z ) = ££ = 2 . ll6pV 2 b 2 sin 8 cos 8 

c. CLT> 


Finally, 


F z = (F z ) x + (P z ) = 3.876pV 2 b 2 sin 8 cos 8 (BIO) 

These results indicate that the fluid within the closed body carries 
45. 4 and the external flow 5k . 6 percent of the z-momentum imparted by 
the airfoil to the fluid. 

It is of interest to compare these results with those of the usual 
momentum treatment of airfoil lift, wherein it is represented that the 

r 

airfoil obtains lift by imparting a downward velocity w = — to a mass 

b 

of air passing through a circle of diameter equal to the wing span 


k2 


L = p — V cos 5 
4 b 

For small values of & (cos 5 -> l), this expression reduces to 

L = 4- pV 2 b 2 sin 5 
o 

so that 

L = 3. 8T6pV 2 b 2 sin 5 

This is identical with the previous result for F z . However, the present 
treatment gives the true physical wake, while the circular form has no 
physical reality. 


i 
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Figure 1. - Trailing vortex system behind a finite- span wing. 


X I V 



Figure 2.- General representation of the vortex- wake inclination. 



Figure 3* - Forces and velocities acting at a section of a finite- span 

wing having high circulation. 
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Figure 6.- Components of the free- stream velocity. 



Figure 7*- Axis system for the initially flat wake strip. 



Figure 8.- Coiling vortex sheet. (See ref. 2.) 
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Figure 9 .- Variation of circulation with radius and the relation of the 

vortex sheet to the core size. 
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Figure 10.- Cross-sectional view of the ultimate wake taken normal to the 

core axes. 
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Figure 11.- Elemental length of a vortex helix of radius r formed by 
the deformation of a vortex sheet wake. 
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Figure 12.- Variation of the function f-]_(r) with radius in the vortex 

core. 
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Figure 13- - General representation of the thickness variation of the vor- 
tex layer behind a finite-span wing. (See ref. 4.) 



Figure lk.~ Addition of the vortex sheet to the vortex core. 



Figure 15.- Force systems for a lifting airfoil. 





Figure 20.- Path of integration for determining kinetic energy external 

to a vortex pair. 


Figure 21.- Geometrical relations for establishing equations ( 64 ) 
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Figure 22.- Geometrical relations for establishing equations ( 67 ) < 
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Figure 2k. - Variation of lift coefficient with aspect ratio and. regions 
in which linear and nonlinear theory apply for calculation of 
induced drag. 




planform and aspect ratio. 






Figure 27.- Geometrical relations for axial- induced velocity field of a 

vortex helix. 



Figure 28.- Ihe streamline i|r = 0 which forms the boundary between the 
open and closed streamline systems. 

NASA - Langley Field, Va. L-1156 


